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! Abstract. Although the cosmic concordance cosmology is quite successfull in fitting data, fine 

i tuning and coincidence problems apparently weaken it. We review several possibilities to ease its 

problems, by considering various kinds of dynamical Dark Energy and possibly its coupling to Dark 



^ |. Matter, trying to set observational limits on Dark Energy state equation and coupling. 
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CN \ 1. Introduction 

^ ■ Until a decade ago two options were in competition: the world could be either SCDM or 

^ ! OCDM. The former cosmology had matter density and deceleration parameters Q-om — 1 

\ and qo = 0.5; the latter had Q.om — 0.2-0.3 and ~ 0. Q-om — 1 was supported by 

■ COBE data and agreed with generic inflationary predictions; — 0.3 was supported 

^ . by evolutionary data and inflation made it acceptable without horrible fine tunings. When 

\Q \ data became too far from SCDM, the mixed model variant [1] became popular. 

O ■ Then SNIa data [2] required qo ~ —0.6-0.7. Q.om could still be ~ 0.3 if the gap up to 

^ . Q.o = I was covered by a substance with p ~ — p, the Dark Energy (DE). The ACDM 

models, considered until then little more than a smart counter-example, begun their fast 
O \ uprise to become the Cosmic Concordance Cosmology (CCC). Then, soon after SNIa 

+l3 ■ data, deep sample analysis [3] and fresh CMB data [4] converged in confirming that 

! Q.om — 0.3 with Q.O — 1 and the CCC became a must. 

\ Only a minority were however happy with A being the cosmological constant. Thus, 

CCC brought the problem of DE nature. Most of this paper deals with DE being a self- 
interacting field, either fully decoupled or decoupled from any other component apart 
5^ \ Dark Matter (DM): dynamical and coupled DE (dDE and cDE), respectively. 

Fig. 1 illustrates why we favor these options. If DE pressure and density meet the 
condition p^e = —pde exactly, then p^e is constant: backward in time, DE rapidly 
becomes negligible and we wonder why just today it became relevant. This is the 
coincidence problem. But, if DE is false vacuum, at the end of the transition, p^/g should 
have became ~ 1 : 10^^ of its pre-transition value. This is a typical ^ne tuning problem. 

The latter problem vanishes in dDE models. In the interaction potential no fine-tuned 
scale appears. However, as shown in panel dDE, the coincidence problem is not eased. 
Then, if DE suitably interacts with DM, densities can evolve as in panel cDE, and also 
the coincidence problem is eased: DE and DM have had comparable densities since long. 

We became accustomed to cosmologies with various components of similar density, 
e.g., baryons and DM. CCC now requires also DE to fall in the same range. This would 
not be an extra requirement if DE and DM are just two different aspects of a Dark Side, 
their interaction being a signature of a common nature. An example is the dual axion 
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FIGURE 1. DM and DE densities in ACDM, dDE amd cDE cosmologies. The solid (dotted) Une shows 
DM (DE) densities in the different classes of models. 



model [5], but here we approach the question from the phaenomenological side. 

We outline soon that DM-DE coupling predicts a significant baryon-DM segregation 
in non-linear structures. Hydrodynamics succeeds in explaining most segregation effects 
observed in the real word, by tuning suitable parameters. If the required tunings conflict 
with other data or other observations require a different behavior of baryons and DM still 
before the onset of hydro, the DM-DE coupling would find an observational support. 

Constraints on cDE from CMB data were first discussed by [6] . In Section 2 the results 
of a fit to WMAPl data [7] are reported. Discussing non-linear constraints is a harder 
task. Even predictions based on the Press & Schechter or Sheth & Tormen [8] formalism 
are hard to obtain. In Section 3 we shall explain why it is so. In Section 4 we shall then 
use ST expressions to predict mass functions. Our conclusions are in Section 5. 

2. Background and linear fluctuations . 

If DE is a scalar field 0, self-interacting through a potential V{^) [9], [10], it is 

Pde = Pk,de + Pp,de = ^^/^a^ + V{^), Pde = Pk,de - Pp,de = W Pde ■ (1) 

Then, if dynamical equations yield pk^del^ ^1/2, it is — l/3^w>— 1. We use 
the background metric ds^ = a^{T){—dT^ + dxidx'); dots indicate differentiation with 
respect to r (conformal time). This DE is dubbed dynamical (dDE) and much work has 
been done on it, also aiming at restricting the range of acceptable w(t)'s, so gaining an 
observational insight onto the physics responsible for the potential V{(j)). 
In this paper we shall consider the potentials 

SUGRA : y (0) = (A«+V0«) exp(4;r0V'«?.) , RP ■ Vi^) = A"+^/(^)" 
([11], [10]), admitting tracker solutions and yielding two opposite w(t) behaviors: 
nearly constant for RP, fastly varying for SUGRA. Then, the effects we find should 
not be related to the shape of V{^) but to the coupling. Most results are shown for 
A = lO^GeV; minor shifts occur when varying logjo(A/GeV) in the 1-4 range. 

DM-DE coupling, fixed by a suitable parameter /3, modifies background equations 
for DE and DM as well as those ruling density perturbations (see, e.g., [6]). 

Linear codes, modified by in this way allow predictions on CMB anisotropics and 
polarization. In this way we fitted dDE and cDE models to WMAPl data [7]. In Figs. 2 
MCMC limits on model parameters are shown, for the dDE and cDE with SUGRA 
potential. The (likelihhod) is 1 .066 (4.7 %) both for ACDM and cDE, 1 .064 (5 %) for 
dDE, with a marginal improvement. More significantly, cDE is not worse than ACDM. 
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FIGURE 2. 1- and 2-a limits on model parameters, obtained from a MCMC fit of dDE and cDE 
(SUGRA) cosmologies with WMAPl. Here X = logjo(A/GeV); A, n^, T are fluctuation amplitude, 
primeval spectral index, opacity. A turns out to be substantially unconstrained. For dDE, is greater 
than for ACDM, but all parameter values are compatible with ACDM values within ^ 2 CJ's. For cDE, the 
coupling j3, at 2-<j, should be <~ 0.3. 



3. Non-linear Newtonian approximation 

While CMB analysis is based on linearized eqs., if we deal with scales well below 
the horizon and non-relativistic particles, the restriction 5c^b ^ 1 can be lifted. An 
alternative Newtonian approximation is then licit, within which tensor gravity and scalar 
field cause overlapping effects adequately described by assuming: 

(i) DM particle masses to vary, so that Mc{x) = Mc{Xi) exp[— C(0 — 0,)]. 

(ii) Gravity between DM particles to be set by G* = jG. 

Here C = A/16;rG/3 /3 , 7 = 1+4/3^/3 (see, e.g., [12]). This approach allowed us to 
study spherical top-hat fluctuations [13], so predicting the cluster mass functions in cDE 
models [14], and to perform n-body simulations [12]. 

The main novel feature of cDE cosmologies we outline in this way is DM-baryon 
segregation, a strong effect, visible in the evolution of a spherical top-hat fluctuation. 

Although in all cosmologies, apart SCDM, the spherical growth must be studied 
numerically, in cDE, just because of the ongoing segregation, the numerical approach 
is essential. In all cases apart cDE, the only variable describing a top-hat is its radius R, 
set to expand initially as scale factor a. Then, the greater density inside the top hat slows 
down R in respect to a, so that the decreasing p{< R) increases in respect to the average 
p. At a time tta, when A = p{< R)/p attains a suitable value x, R starts decreasing, to 
formally vanish within a finite time tc- For sCDM, x = (3;r/4)^ and tc/tta = 2. In other 
non-cDE models, x and tc/tta take slightly different values. 

However, when p(< R) increases, unless the heat produced by the pdV work is radi- 
ated, virial equilibrium is soon attained, and any realistic fluctuation stops contracting at 
a radius Ry. In sCDM, Ry/Rta = 1/2 and, taking into account the symultaneous growth 
of a, the virial density contrast is then Ay = 32;^ ~ 180, indipendently of tc- Values are 
slightly different for other ACDM and dDE and a usual assumption is that the system 
relaxes into its virialized configuration just within tc. 
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FIGURE 3. Density profiles at different a values for j3 = 0.2 model. Notice the progressive deformation 
of the baryon profile (dotted lines) in respect of the DM profile (solid line). 



Av values for ACDM and dDE were obtained by [15] and [16, 17], respectively. The 
same problem was treated in [13] for cDE. In this case, and Rc (baryon and DM 
radii) initially grow as the scale factor, but Rj, soon exceeds Rc, because of the stronger 
effective gravity; peripheral baryons then leak out from Rc, so that DM no longer feels 
their gravity, while baryons above Rc feel the gravity also of external DM layers. Then, 
above Rc, the baryon profile is no longer top-hat, while a fresh perturbation in DM 
arises. Re-contraction will then start at different times for different components and 
layers. Similarly, virialization conditions are fulfilled earlier by inner layers, although 
outer layer fall-out shall later perturb them so that the onset of virial equilibrium is a 
multi-stage process. Furthermore, when the external baryons fall-out onto the virialized 
core, richer of DM, they are accompanied by DM materials originally outside the top- 
hat, perturbed by baryon over-expansion. 

Each layer and substance, in such system, feels a different force; each shell needs then 
to be considered separately; this is why the numerical problem is far more intricate. 

The collapse stops at virialization. Fig. 3 assumes the growth to stop when all DM 
originally in the top-hat, and baryons inside it, virialize. Already at a = 0.2, before 
the turn-around, the baryon top-hat boundary is no longer vertical. At a = 0.6 (~turn- 
around), a similar effect is visible also for DM. The effect is even more pronounced at 
avir — 0.92. The same effects, just slightly weeker, are present also for j8 = 0.05 

Fig. 3 shows also that baryons leak out from the DM top-hat. For a ~ 0.92, when 
DM and inner baryons virialize, ~ 10% ( ~ 40%) of baryons are out of the top-hat, if 
j8 = 0.05 (j8 = 0.2). These values increase by an additional 10% in the RP case. 

The virialization of materials within R requires that 2T{< R) = RdU{< R)/dR. 
Kinetic and potential energies have fairly straightforward expressions. The only point 
to outline is that DM-DE energy exchanges, for the background, are accounted for by 
background evolution eqs., so that, when fluctuations are considered, the background 
contribution must not considered again. 
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4. Mass functions in cDE theories 

Let us compare the actual growth of a top-hat fluctuation with its growth if we assume 
linear equations to hold, indipendently of its amplitude. While the real fluctuation 
abandons the linear regime, tums-around, recontracts and recoUapses at the time T^o 
linear equations let that fluctuation steadily grow, up to an amplitude 5rc at the time Xrc- 
The linear evolution does not affect amplitude distributions. If they are initially dis- 
tributed in a Gaussian way, at Xrc we can still integrate the Gaussian from 5rc, so finding 
the probability that an object forms and virializes. This is the basic pattern of the PS-ST 
approach, that we apply also to cDE cosmologies. For them, however, Xrc is different for 
DM and baryons. Viceversa, if we require To to coincide with Xrc, we have two different 

linear amplitudes, 5rc and 5rc , yielding recoUapse at t„ for all baryons or all DM. 

The j8 dependence of 5^^ at z = and their z dependence are given in [14]. Setting 
^m — ^m/<^m, the PS differential mass function then reads 

YiM) = y/2/jt{p/M) / dVM (dVM/dM) VMexp{-v^/2} , (2) 

J 8rcl Om 

using here 5rc\ 5rc \ or any intermediate value, according to the observable to be fitted. 
The ST expression is obtainable from eq. (2) through the replacement 

VMexp(-f )^^'v'M(l + v'^^/^)exp(-^), with ^' = 0.322, v'J, = 0.707 v^, 
meant to take into account the effects of non-sphericity in the halo growth. 

In ACDM or dDE, the mass M in ST expressions is the mass originally in the top- 
hat. In cDE, indipendently of the value taken for 5cr, virialized systems will be baryon 
depleted. In fact, mass function built by taking 5^^, concern objects before a part of the 
initial baryon content has fallen out. But, even if we wait for a total or partial baryon fall 
out, by using 5^^ or an intermediate value, the initial baryons come back carrying with 
them DM layers initially external to the top-hat. A prediction of cDE theories, therefore, 
is that Q-c/O-b, measured in any virialized structure, exceeds the background ratio. 

We shall plot mass functions obtained using either 5c^rc or ^c,rc- Actual data should 
fall in the interval between these functions. If, during the process of cluster formation, 

outer layers were stripped by close encounters, data shall be closer to the 5c'^)c curve. 

Figure 4 shows n{>M) — dM' \]/{M') obtained through eq. (2). The lower panel 
shows the ratio between halo numbers for each model and ACDM. No large differences 
between ACDM and dDE are found at z = . Shifts are greater between ACDM and cDE. 
For j8 = 0.20, clusters with M >~ 3 ■ lO^^/j-iM© are half than in ACDM. For /3 = 0.05 
the shift is smaller, hardly reaching 20%, still opposite to dDE. The r.h.s. panel shows 
that these features are true also for RP. More RP results can be found in [14]. 

To discriminate between models, cluster evolution predictions should provide num- 
bers per solid angle and redshift interval (rather than for comoving volumes, see [18]) as 
in the upper panels of Fig. 5. In the lower panels, SUGRA models and ACDM are com- 
pared. We consider two masses: 10^'^ and 4 • lO^^h^^MQ-, for the latter scale, a magnified 
box shows the expected low-z behavior showing how we pass from numbers smaller 
than ACDM to greater numbers, for /3 = 0.05, at a redshift z — 0.7. The high-z behav- 
iors for j8 = 0.05 or 0.2 lay on the opposite sides of ACDM. 
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FIGURE 4. Cluster mass function at z = 0. Four models are considered: ACDM, SUGRA dDE with 
A = lOOMeV, and two SUGRA cDE with J3 = 0.2 and 0.05. Dashed areas are limited by the mass 

(c) (b) 

functions for 5c,rc or 5c,rc in cDE models. The lower panel yields ratios in respect to ACDM. The r.h.s. 
panel overlaps cDE results for SUGRA and RP. 




FIGURE 5. Number of clusters with M > lOi'^ft-'M© (left panel) or 4 • lO^'^h'^MQ (right panel) in a 
fixed soUd angle and redshift interval. The upper linear (lower log) plot outiines shifts at small (large) z. 



5. Conclusions 

In this paper we discussed individual cluster formation and cluster mass functions. A first 
finding is the expected baryon-DM segregation causing baryon depletion in clusters. 
Depletion could be even stronger if the outer layers are stripped out by close encounters 
during the formation process. Preliminary results of simulations confirm these outputs. 

As far as mass functions are concerned, a first significant feature is that the discrep- 
ancy of dDE from ACDM is partially or totally erased by a fairly small DM-DE cou- 
pling, and many cDE predictions lay on the opposite side of ACDM, in respect to dDE. 

At z = 0, a shortage of large clusters is expected in cDE. Therefore, if ACDM is used 
to fit data in a cDE world, cluster data yield a smaller than galaxy data. 
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When we consider the z dependence, we see that (i) when passing from ACDM to 
uncoupled SUGRA, the cluster number is expected to be greater, (ii) When coupling is 
added, the cluster number excess is reduced and the ACDM behavior is reapproached. 
(iii) A coupling j8 =0.05 may still yield result on the upper side of ACDM, while 
j8 = 0.2 displaces the expected behavior well below ACDM. 

This can be fairly easily understood. When p^e keeps constant, while Pm°^ (1 
DE relevance rapidly fades. In dDE models, instead, p^e (slightly) increases with z and, 
to get the same amount of clusters at z = 0, they must be there since earlier. Coupling 
acts in the opposite way, as gravity is boosted by the field. In the newtonian language, 
this means greater gravity constant (G*) and DM particle masses at high z, speeding up 
cluster formation: a greater /3 needs less clusters at high z to meet their present number. 

The overall result we wish to outline, however, is that non linearity apparently boosts 
the impact of coupling so that there are quite a few effects from which the coupling be- 
tween DM and DE can be gauged. Most of them are just below the present observational 
threshold and even slight improvements of data precision will begin to allow discrim- 
inatory measurements. An essential step to fully exploit such data will be performing 
n-body and hydro simulations of cDE models. Through this pattern we can therefore 
expect soon more information on the actual nature of the Dark Side. 
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